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1. Introduction 

We started the session by using elimination to convert a first order linear 
system to an ordinary differential equation in one of the variables. In this 
note we'll see how to go the other way. That is, to convert a second order 
ODE to a 2 X 2 system of first order equations. Looking at a second order 
system this way gives us an important way to visualize the second order 
equation. 

Since elimination takes us in the other direction we might call this pro- 
cess anti-elimination. This term is not entirely standard, but it will serve 
us nicely. 

We won't cover it in this course, but every method we saw for numerical 
solutions to first order ODE's goes through without change to systems of 
first order equations. Thus, anti-elimination allows us to use numerical 
techniques on second order ODE's. Indeed it can be used to convert any 
order ODE to a system of first order equations. 

2. Anti-elimination and the Companion Matrix 

Suppose we have a second order homogeneous linear equation, say 



We can derive a first order linear system from this, by using the following 
trick. Introduce a second variable defined by 



X + bx + kx = 0. 



(1) 



y = x 



Substituting y = x and y = x into equation (1) we get 

y + by -\-kx = Q ^ y = —kx — by. 



We now have a first order system 



X 



y 



V = 



kx 



- by 



The corresponding coefficient matrix is 
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This is called the companion matrix of the equation (1). In this case, the 
solution vector u{t) = ( ^^^^ J . It records both the solution to (1) and its 
derivative. 

Example 1. Consider the equation x — x + = 0. The companion matrix 
1 



IS 



-5/4 1 



What do solutions of this system look like? The characteristic polyno- 
mial of the second order equation is 

p{s) = - s + 5/4 = (s - (l/2))2 + 1. 

So, the roots are r = (1/2) ±i. From unit 2, the general solution in amplitude- 
phase form is given by 

x{t) = Ce^^^cos{t -(p), 

where C and cp are constants. These oscillate under an exponentially grow- 
ing envelope. The derivative does the same, but is off phase. This means 
that the trajectory traced out by {x,y) = (x,x) is an expanding spiral. 

Example 2. (Elimination followed by anti-elimination) 
Earlier in this session, we learned how to solve systems by elimination. 
What happens when we do elimination followed by anti-elimination? Let 
us re-visit the example about Romeo and Juliet. The second order system 
describing their mutual feelings was 

R'= \r + J (2) 

J' = -Ir^Ij (3) 

Let us first eliminate /, to get a second order equation in R. From (2), 
J = R' -{1 /4)R. Substituting this into (3) gives 

R"-Ir' = -^R+^R'-^R ^ R"-R'-^R = 0. (4) 
4 16 4 16 4 

(We saw this equation in example 1.) Applying anti-elimination gives 

R' = Y 

y = iR + y. 
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This is a different system from the one we started with. The companion 
matrix of the ODE (4) is different from the original matrbc associated to the 
system. 
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